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Abstract 

The Renormalization Group equation describing the evolution of the metric 
of the nonlinear sigma model poses some nice mathematical problems involv- 
ing functional analysis, differential geometry and numerical analysis. In this 
article we briefly report some results obtained from the numerical study of 
the solutions in the case of a two dimensional target space (deformation of the 
0(3) cr-model). In particular, our analysis shows that the so-called sausages 
define an attracting manifold in the [/(l)-symmetric case, at one loop level. 
Moreover, data from two-loop evolution are used to test the association 
between the so-called SSM V field theory and a certain £/(l)-symmetric, fac- 
torized scattering theory (FST). 
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I. RENORMALIZATION GROUP EQUATION 



The perturbative renormalization of the non linear a-model 
of the metric according to the (one-loop) equation 



dg 



dt 



■—R lJ+ 0(R 2 ). 



gives rise to a deformation 



(1) 



This second order non linear partial differential equation (PDE) has been studied in the 
simplest case of 2-dimensional target manifolds in ref. Jl|. A whole family of solutions is 
known for the topology of the sphere S 2 or for the torus. In this letter we consider the case 
of S 2 only, but our method is easily adapted to the toroidal case. As is well known we can 
introduce local coordinates {y, p} in such a way that the metric is conformally flat: 



9ij ^ Oij 

and the RG equation valid to all loop reduces to a single non linear PDE, 



(2) 



where (3{R) is the exact (3— function and 

R = -e 



d 2 d) <9 2 < 



d 2 y d 2 (p 



(3) 



(4) 



is the scalar curvature. 

A family of solutions in the case of one-loop, namely /3(x) 
in These can be obtained starting from the ansatz 



-x/47r has been presented 



log[a(t)+b(t)f(y) + c(t)h(<p)} . 



(5) 



The only solutions of this kind which can be extended to t — > — oo (UV limit) without 
encountering singularities and which exibits a residual U(l) symmetry are the so-called 
sausage solutions, parametrized by a single real constant v: 



<f>{y,t) = -log{a{t) + b{t) cosh 2?/) 



(6) 
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These trajectories in the space of all metrics are believed to be the one-loop approximation 
of some integrable deformation of the costant-curvature trajectory (0(3) a-model), to which 
they tend as v — > 0. To discuss the general solution of the RG equation we have to construct 
a numerical integration algorithm. The first problem is posed by the divergence of the 
conformal factor as \y\ — ► oo: <p ~ — 2\y\ as \y\ — > oo. This general property (which 
follows directly from Gauss-Bonnet theorem) causes the factor exp(— <p) to diverge, hence 
amplifying the numerical error of the differential term at large \y\. To overcome this difficulty 
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we introduce a background field (po(y,t), and we consider the equation for the shifted field 
= 7] + 0o- The background field is conveniently chosen as the constant curvature solution 



log 



A(t) 



An cosh 2 y 



(7) 



where according to the RG equation A(t) satisfies 

A(t)A(t)~ l = ^TtAit)- 1 ) . (8) 
Introducing a new evolution parameter ("time") r, defined by 

A(t) = Ane- T (9) 
the RG equation expressed in terms of the shifted field rj reads 

with 

R( V ) = R(<f> + V ) = e T -"(2 - A oV ) . (11) 

Here A = cosh 2 ydy + 9 2 is the standard 0(3)— invariant Laplacean on the sphere. The 
original time scale can be simply recovered by quadratures 

rr dr' 

( ° - * = L wf) ■ (12) 

At this stage we restrict our attention to U(l)— symmetric solutions r\ = rj(y,t). 

II. NUMERICAL INTEGRATION OF THE RG EQUATION 

At this point it is convenient to adopt the standard spherical coordinates by letting 
y = log(cot($/2)). To get a good accuracy in the evaluation of the Laplacean we apply 
the spectral method, that is we expand r\ in Legendre polynomials P n (cos?9), which are the 
eigenfunctions of A . Thus, we developed the following finite-dimensional implementation 
of the standard expansion in Legendre polynomials. 

Let be the zeroes of the L-th Legendre polynomial Pl{x); we adopt {x^ \ j = 

1, . . . , L) as our finite gridQ and sample the field rj(y, t) at the image points y^ L \ 
The finite Legendre expansion is then realized as follows: 



Solutions with reflection simmetry (x — > —x), such as the sausage solutions, can be studied by 
a restriction to even-order polynomials; our algorithm can be easily adapted to this case, with a 
gain of a factor four in memory requirements 
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oo L— 1 

v( x ) = ^2viPi( x ) tt^rjiPiix) (13) 
;=o 1=0 

Vi = \ J\ dx r ? (x)P z (x)(2Z + 1) « i Vf >(2Z + 1) = (^)j (14) 

where are the Gaussian integration weights || for Legendre polynomials. This finite 
expansion allows us to represent the Laplacean exactly on polynomials of degree less than L, 
since the coefficients r\\ are exact in this case. The lack of a fast implementation analogous 
to FFT limits our algorithm in practice to L ~ 200 on current workstations, but this proves 
to be adequate for our purposes. 

Having determined the finite transform £ on the basis of Legendre polynomials, the 
action of the Laplacean A is represented by a matrix £~°°*£, where A is diagonal, with 
eigenvalues {—1(1 + 1) | I = 0, . . . , L — 1}. In terms of this representation is quite easy to 
compute the spectrum of zero-modes of the field 0, a problem considered in ref. [|IJ. We 
just have to diagonalize the finite matrix — (1/2) A + (1/8)R, where A = exp(r — 77) A . The 
results agree with the previously computed ones for the sausages J7J (notice however that 
the present method is much simpler and the spectrum can be computed in parallel with the 
RG evolution). 

With the finite Legendre-transform algorithm at hand, we can now consider the integra- 
tion of eq. (|i~0"l). We can work to any loop, provided we know the corresponding (3— function. 
We have implemented the algorithm in matlab || which provides efficient routines of dia - 
go -nalization and of adaptive-step integration. The accuracy of the code has been tested 
on the known one-loop sausage solution, attaining a typical maximal deviation of 1 part 
in 10 9 over a time interval —3 < r < 2 and v < .25. In Fig.l we show a typical one-loop 
sausage evolution; to aid the visualization we have recovered the embedding of the surface in 
R 3 in such a way that the induced metric coincides with exp(</>)5y. The accuracy is limited 
essentially from the large eigenvalues of the Laplacean which grow as the square of the finite 
grid dimension. Moreover as v and/or — r increase the curvature tends to be confined at the 
extremities of the sausage, which requires finer and finer grid. Presently, up to 200 points, 
we cannot go beyond r « —4 for v ~ .25, but there is no limit in principle. 

The algorithm can now be applied to investigate, at one-loop, the existence of attracting 
manifolds, in the space of all metrics. A conjecture of Fateev and Zamolodchikov states that 
the sausages constitute a stable manifold and all other geometries converge to some sausage, 
parametrized by the real v, in the IR direction r — > 00. 

This fact manifests itself quite clearly in our numerical data. We may define a distance 
function by 

dist(7 7l , 7 72 ) = J dx[R(7 7l ) - R(t? 2 )] 2 (15) 

and measure the distance to the sausage hypersurface by 

D(t) = Mdist( V (T), Vv (r)) (16) 

where r\ v is the sausage solution specified by v. The data show a clear exponential decay 
D(r) ~ ylexp(— mr), with m ~ 8 regardless of the initial surface (see Tab.l). For any given 
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starting geometry we record the value v e g where the infimum is reached; its limit iA; m as 
r — > oo gives a definition of sausageness of any given surface. For instance, given an ellipsoid 
with cylindrical symmetry and eccentricity e we can measure e(h>\ im ), at one-loop order (see 
Fig. 2). This shows a remarkable property: e is a universal function of s = vu m ^j-, where f 
is the starting scale of the evolution. This function seems to be e = tanh s, at least for small 
s. 



III. VARIATIONAL EQUATIONS AND STABILITY 

Another way to discuss the attracting nature of the sausage manifold is to study the 
Jacobi variational equations around the sausage solution. The linearized equations take on 
the following form ( calling x the variation): 

d rX = H( V ,r) X (17) 

where 

H(v,r) = ^0y{ eT - VA o + R(v)) • (18) 

The spectrum of H is not a priori of much significance, since the evolution equations are 
time-dependent. However, if we rely on the adiabatic approximation, the spectrum is di- 
rectly related to stability. Applying the finite Legendre-transform £ the spectrum of H 
is easily reduced to that of a Hermitean matrix; choosing for 77 some sausage r\ v we find a 
single positive eigenvalue denoting an unstable mode. As we presently show, this is easily 
interpreted. 

It is immediate to find that the Jacobi variational equation at one loop: 

drX = \ (exp(-^)Ao + e- T RM) X (19) 

admits the two solutions 

Xi = R(Vu) (20) 
X 2 = 1 - (1/2) eM~r)RM- (21) 

Their existence is not surprising: the family of sausage solutions ([)]) depends on two param- 
eters t and v, so that the derivatives with respect to t and v give rise to two independent 
solutions of the variational equations. Now Xi, which comes from the time derivative rep- 
resents the unstable solution; this kind of instability is of no concern, since it corresponds 
to a simple redefinition of the initial time parameter t and it can be fixed by restricting to 
a given initial area. The other solution \2 comes from the ^—derivative and is tangent to 
the manifold of sausages. The component of a generic variation x which is orthogonal to 
both Xi an d X2 provides a measure of the distance of a generic solution from the sausage 
manifold. 
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Numerical results (Tab. 2 and 3) show that the orthogonal part is exponentially decreasing 
with a slope (~ 8.73)0 much higher than the slope of Xii which from equation (^l|) is found 
to be 2. The point is that there is an overall convergence in the infrared toward the constant 
curvature metric. All sausages converge to a sphere with vanishing radius, forcing the decay 
of the mode xi\ however, the rate of this convergence is slower than the rate of decay of all 
other modes. Fig. 3 shows how the longitudinal component (plotted in the vertical direction) 
is stil large when the transversal one is already negligible. In Fig. 4 there is an impressive, 
though only qualitative proof of what we found (one should compare Fig. 4 with Fig.l, where 
the sausage evolution is plotted). 

Conversely, the evolution in the ultraviolet direction r — > — oo is strongly unstable, which 
makes it quite hard to follow numerically; all truncation errors are chaotically amplified and 
the calculation becomes rapidly unreliable. 



IV. THE FREE ENERGY OF THE SAUSAGE SIGMA MODEL 

In ref. [TJ the model corresponding to U(l) — symmetric, one-loop deformation of the 0(3) 
sigma model (termed SSM U model) is associated with a U(l) — symmetric deformation of the 
0(3) FST. To test the correctness of the identification, some physical quantities calculated 
from S"-matrix data via Thermodynamic Bethe Ansatz tecniques are compared in the UV 
region with the corresponding quantities obtained perturbatively from the one-loop sausage 
action 



A 



SSMv 



J 1 3 H + LcoshGf ' 



Their non-trivial matching is a strong (though indirect) proof of the correctness of the 
association. 

However, while for the FST we have exact (at least in principle) results, for the field 
theory we are strongly limited by the fact that the "sausage metric" is a solution to the RG 
equation only to one-loop. So, it would be very important to have higher-loop solutions of 
the RG equation with the same features of the sausage, so to confirm (or deny) the validity 
of the association proposed in |TJ. This is not to be taken for granted. To any relativistic 
FST one would like to associate an integrable QFT; however, even if one assumes that the 
action ( |2"2| ) is the one-loop approximation of an exact trajectory of the RG group, which 
would define the exact sausage model, there is no a prori reason to expect such model to 
be integrable. In particular, S 2 with any {/(l)-symmetric metric is not a symmetric space 
and one cannot use the results of refs. [|J to establish the classical and/or quantum 
integrability of the model. 

Considering the zero-temperature free energy f(A) in the presence of a constant external 
field A, the one-loop calculation based on the action ([22] ) gives the result JlJ 



2 This slope is somewhat higher than the one found in sect. II; there we were actually considering 
variations with a nonzero slowly-decaying component along \2- 
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f(A) = -A 2 e^ ^ , (23) 

where to — t = log A/ A, A is a subtraction scale and one works in the "scaling limit": 

v — > , t — > — oo , z4 fixed (24) 

This limit serves to eliminate all higher-loops contributions, which cannot be properly taken 
into account using the one-loop action (|2"2"|). To obtain higher-loop correction to ( |2"3"D one 
needs the corresponding higher-loop solution of the RG equation. Formula (|23|) however 
remains valid, assuming as natural that the higher-loop action is an even, convex function 
of A. 

With our algorithm, the quantity in (|23|) is evaluated as 

f(A) num = —A 2 e 7l( - t ' Xinad ^e' t ' o( - t ' y=0 ^ = A 2 e n( ' t ' Xined ^ — — - (25) 

4:71 

where x me d is the middle point of the finite grid with represents the coordinate x = cos 9. 
The two-loop /3— function is well known 0: j3(x) = —x/4n — x 2 /(4tt) 2 . Lacking an exact, 
£7(1) — symmetric two-loop solution of the RG equation, the first thing one could try is a 
linerization around the one-loop sausage solution, regarding v as perturbative parameter. 
Notice in fact that v may be scaled out from the sausage solution, through the substitutions 

0^0 + log£-, t^vt (26) 

which leave the one-loop equation invariant (actually, with our choice, the scalar curvature 
R is conveniently scaled by a factor of 4tt). On the other hand the two-loop equation is 
transformed by Q2"6D to 

& f t =-R- vR 2 (27) 

where the last term could now be treated perturbatively for small sausage deformations. 
However, even the linearized equation proves to be too difficult to solve exactly, when the 
correct boundary conditions for <p are taken into account, and we resort to a numerical 
approach. 

Our idea is the following: in the scaling limit (far UV and small deformations), the exact 
solution must have the same features of its various loop-wise approximations. We can now 
study with our algorithm the two-loop evolution, using as initial data in the far UV (say at 
t = t) a metric with the "sausage features " (e.g. ge itself, one-loop solution). We 

then evaluate our physical quantity in a fixed scale interval (^1,^2) and we compare it with 
the same quantity evaluated starting from the FST by means of Bethe Ansatz techniques 
||. We expect that the difference between the two quantities tends to zero when the starting 
scale of the evolution is pushed farther and farther in the UV. At one-loop level this is easily 
verified, using as starting data, e.g., the function 

0(t, y) = 4> (t, y) + log - - log(l + 2e u{I - to)/27T cosh 2y + e ^- to)/7T ) (28) 
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The above mentioned difference vanish exponentially when the initial scale t is sent to — oo. 

In ref. Q the Bethe Ansatz techniques are used to compute the free energy as a function 
of the chemical potential h coupled to the physical particles. The input is the factorized, 
U(l)— symmetric S-matrix with deformation parameter A relative to the S"-matrix of the 
0(3) sigma model. The output reads, in the scaling limit A — > 0, h — ► oo, Alog/i fixed, 



m = -fl^ 

Z7rA 1 + U 

2A/(1-A) 



l + 4A-^logi-^ + 0(A 2 log 2 A) 
1 — u A IX 



(29) 



where u = ( me ^ ) , and m is the mass of the physical particles. The connection 

among the various quantities in (|29"D and those contained in the perturbative expressions is 
established as follows: 

TL- = A + 0{u) (30) 
h = A + 0{v) (31) 
± = ±.[l + 0{u)] (32) 

The relation between the "scale variables" z = v(to — t)/2n and u follows by noticing that 
in the perturbative calculations we put t — t — log(A/A). Substituting ( [30|) -( j32]) in (|29"D we 
obtain 



A 2 " 7T 



- tanh(,/2) + - J^_ - M2*q - e-«)) 
i/ cosh 2 (^/2) cosh 2 (z/2) 



(33) 



where the subscript "an" stands for "analytical" (compare (p5|)). The first term comes 
from the one-loop action (|22]), while the second and the third represent the two-loop 
contributions. The presence of the non-analytic behaviour ulogu ensures the correct 
0(3)— symmetric limit {y — > at fixed t — t). Corrections of order v are not known; 
to establish them we should know the corrections of the same order to eqs. (|30D - (132]) . Our 
aim is now to compare expressions f(A) aii and f(A) num . 

We have collected numerical data relative to the two-loop evolution of sausages with v 
fixed, starting at different initial scales f (recall that r is related to to — t by fljjD). Fig. 5 
reports the results; here v = 0.1 and the r interval is chosen as (r x = — 3.85,r 2 = —3.8). 
The quantity A is defined as 

A(f)= max \f(A) & JA 2 -f(AUJA 2 \ (34) 

T1<T<T2 

and it is plotted for different values of the initial scale f. 

A does not decrease when the scale is pushed towards far UV; instead, initially increase and 
then reach an asymptotical value. The value to which A tends is of order 10~ 2 ; that means 
a discrepancy of about 5 parts in 10000 relative to the value of the free energy. This is well 
above the numerical errors of our algorithm. 

We controlled if this discrepancy can be reabsorbed in a "renormalization" of u; leaving it 
as a free parameter in (B3|), we seeked for a minimum of A in a neighborhood of 0.1. The 
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result was negative. 

In principle, the non-zero asymptotic value of A may be eliminated by the 0{v) correction; 
to test this hypotesis we repeated the evaluation of A for different v values, and corre- 
spondingly redefined the r— interval so to leave the interval (uti, vt-i) fixed, according to the 
scaling limit. The results are summarized in Fig. 6: the asymptotic value of A seems not 
to depend appreciably upon v. It remains always of the same order of magnitude, and does 
not grow with v. 

V. CONCLUSIONS 

The RG equation for the non linear sigma model with fields taking values on a target 
space with the topology of the 2-sphere is a nonlinear partial-differential equation, possi- 
bly nonlocal (if the (3 function were known to all loop). These features makes it hard to 
caracterize the general solution in an analytical way. The numerical approach appears at 
the moment to be the only way to obtain quantitative informations about the higher-loop 
solutions. 

The spectral method presented in this letter proved to be a useful tool for the construction 
of an efficent algorithm; this remains true even if the exact RG equation should turn out 
to be nonlocal. Another future developement concerns the study of solutions without any 
residual symmetry, hopefully with the help of some "fast" algorithm analogous to the well 
known FFT for the Fourier case. 

In the one-loop approximation we have shown the attractive nature of the 
U(l)— symmetric family of solutions, the so-called sausages. 

As far as the two-loop evolution is concerned, however, our data show a discrepancy 
between calculations based on the RG equation and analogous calculations based on the 
U(l)— symmetric FST proposed in [[j]. 

This probably means that the identification put forward in [[]]] is not correct, in the sense 
that the U(l)— symmetric FST after all does not correspond to a U(l) — symmetric nonlinear 
sigma model reducing to the sausage at one-loop. On this point, the following alternative 
conjecture seems natural. 

The factorized S— matrix put forward in |I[ is a quantum-group deformation of the 
0(3) S— matrix, in the sense of, e.g., ref. ||. In practice, the triplet of massive particles 
of the model forms a spin-one irreducible representation of SU(2) q and scatter through a 
SU (2) q — invariant S— matrix with q ~ e tX . Hence, the QFT corresponding to such FST 
should enjoy a nonlocal hidden SU(2) q invariance, whose U(l) subgroup coincides with the 
manifest U(l)— symmetry locally implemented. In the usual spirit of the sigma models, in 
which the field-theoretical symmetries follow from geometrical properties of the target mani- 
fold, we expect this quantum-group invariance to follow from the noncommutative geometry 
of a q— deformed target manifold SU(2) g /U(l) (attempts in constructing q— deformed sigma 
models can be found in refs. [Kj fllTf ). In other words, we suspect that the FST put forward 
in does not correspond, beyond the one-loop approximation, to any conventional nonlin- 
ear sigma model (that is to a model in which the fields take values in an ordinary manifold 
with commuting geometry). It would be interesting to explain why there exist a one-loop 
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conventional sigma model, perhaps through a suitable expansion near q — 1. This matter 
clearly requires further investigations. 
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FIGURES 




v = 0.2 




FIG. 3. Three-dimensional view of the decay of the longitudinal and transversal modes 




FIG. 4. A pictorial way of illustrating the attractive nature of the sausage solution 
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FIG. 5. The difference, for various values of r, between the free energy evaluated from RG 
equation and the corresponding quantity evaluated from FST data 
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FIG. 6. The asymptotic value of A for different values of the parameter v 
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TABLES 



TABLE I. Evolution of an ellipsoid; D is the distance to the sausage manifold. 



r 




D 


-4.00000 


0.03761 


0.04265 


-3.90000 


0.03609 


0.01453 


-3.80000 


0.03557 


0.00551 


-3.70000 


0.03536 


0.00217 


-3.60000 


0.03528 


0.00088 


-3.50000 


0.03524 


0.00056 


-3.40000 


0.03523 


0.00015 


-3 30000 


03599 


00006 


-3.20000 


0.03522 


0.00003 


-3.10000 


0.03522 


0.00001 


-3.00000 


0.03521 


0.00000 


-2.90000 


0.03521 


0.00000 


-2.80000 


0.03521 


0.00000 


-2.70000 


0.03521 


0.00000 


-2.60000 


0.03521 


0.00000 


-2.50000 


0.03521 


0.00000 
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TABLE II. Decay of longitudinal component (x\\) and of the transversal one (x±) at v = .1 



r 


X\\ 


x± 


-1.9000 


0.0047 


0.1772 


-1.8000 


0.0055 


0.0734 


-1.7000 


0.0051 


0.0303 


-1.6000 


0.0043 


0.0125 


-1.5000 


0.0036 


0.0051 


-1.4000 


0.0030 


0.0021 


-1.3000 


0.0025 


0.0009 


-1.2000 


0.0020 


0.0004 


-1.1000 


0.0017 


0.0001 


-1.0000 


0.0014 


0.0001 


-0.9000 


0.0011 


0.0000 


-0.8000 


0.0009 


0.0000 


-0.7000 


0.0008 


0.0000 


-0.6000 


0.0006 


0.0000 


-0.5000 


0.0005 


0.0000 



TABLE III. Decay of longitudinal component (x\\) an d of the transversal one (x±) at v = .3 



r 


XII 


X± 


-1.9000 


0.0761 


0.4415 


-1.8000 


0.0963 


0.2057 


-1.7000 


0.0944 


0.0939 


-1.6000 


0.0848 


0.0421 


-1.5000 


0.0734 


0.0186 


-1.4000 


0.0623 


0.0081 


-1.3000 


0.0524 


0.0035 


-1.2000 


0.0438 


0.0015 


-1.1000 


0.0364 


0.0006 


-1.0000 


0.0302 


0.0003 


-0.9000 


0.0250 


0.0001 


-0.8000 


0.0207 


0.0000 


-0.7000 


0.0171 


0.0000 


-0.6000 


0.0140 


0.0000 


-0.5000 


0.0116 


0.0000 
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